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Abstract 



The origin of classical predictability is investigated for the one dimen- 
sional harmonic chain considered as a closed quantum mechanical system. 
By comparing the properties of a family of coarse-grained descriptions of the 
chain, we conclude that local coarse-grainings in this family are more useful 
for prediction than nonlocal ones. A quantum mechanical system exhibits 
classical behavior when the probability is high for histories having the cor- 
relations in time implied by classical deterministic laws. But approximate 
classical determinism holds only for certain coarse-grainings and then only if 
the initial state of the system is suitably restricted. Coarse-grainings by the 
values of the hydrodynamic variables (integrals over suitable volumes of den- 
sities of approximately conserved quantities) define the histories usually used 
in classical physics. But what distinguishes this coarse- graining from others? 
This paper approaches this question by analyzing a family of coarse-grainings 
for the linear harmonic chain. At one extreme in the family the chain is di- 
vided into local groups of N atoms. At the other extreme the N atoms are 
distributed nonlocally over the whole chain. Each coarse-graining follows the 
average (center of mass) positions of the groups and ignores the "internal" 
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coordinates within each group, these constituting a different environment for 
each coarse-graining. For an initial condition where long wavelength modes 
are excited and short wavelength modes are distributed thermally we find that 
the coarse-grained positions obey deterministic equations of motion accom- 
panied by noise. The noise is greater the more nonlocal the coarse-graining. 
Further, the deterministic equations require more time steps to evolve over a 
given time interval for the nonlocal coarse-grainings than for the local ones. 
A continuum limit is possible only for the near local coarse-grainings. For 
parameters of the model characteristic of realistic situations these features 
strongly favor the local coarse-grainings over the nonlocal ones for prediction. 
Each of these differences can be traced to the approximate conservation of 
the local center of mass momentum. We then consider the chain quantum 
mechanically and show that for realistic parameters, all the coarse-grainings 
decohere rapidly compared to dynamical time scales. We conclude that noise, 
decoherence, and computational complexity favor locality over nonlocality for 
deterministic predictability. 
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I. INTRODUCTION 



As far as we know them, the laws of physics that apply universally to all physical systems 
are quantum mechanical. The universe at a fundamental level is therefore characterized 
by indeterminacy and distributed probabilities. The wide range of applicability of classical 
deterministic laws is an empirical fact to be explained from the universe's quantum dynamics 
and initial quantum state. This paper investigates the origin of classical predictability 
for the very simple model of a linearized chain of idealized atoms in the context of the 
quantum mechanics of closed systems fl|-|3| , most generally quantum cosmology. We exhibit 
decoherent sets of coarse-grained histories for which, given suitable restrictions on the state, 
the probability is high for histories exhibiting the correlations in time governed by the 
classical wave equation. We shall compare these quasiclassical coarse-grainings with a class 
of others and analyze why the classical coarse-grainings are the most predictable in the class. 

Why do we raise the question of the origin of classical predictability over seventy years 
after the initial formulation of quantum mechanics? Every quantum mechanics text contains 
some treatment of this question. Ehrenfest's theorem is the starting point for one such 
discussion. For a nonrelativistic particle of mass m moving in one dimension in a potential 
V(x), Ehrenfest's theorem is the exact relation between expected values: 

This is not a deterministic equation of motion, but for certain states, typically narrow wave 
packets, the expected value of the force may approximated by the force at the expected 
value of the position, thereby giving a classical equation of motion for that expected value: 

d*(x) dV((x)) 
m dt* ~ dx ■ [ ] 

This equation shows that the orbit of the center of a narrow wave packet obeys Newton's 
laws. 

This kind of elementary derivation already exhibits two necessary requirements for a 
quantum system to exhibit classical deterministic behavior. Some coarseness in the descrip- 
tion is needed, as well as a restriction on the initial quantum state. However, otherwise this 
kind of demonstration does not address the issues we hope to discuss in quantum cosmology 
for the following reasons: 

• The behavior of expected values in time is not enough to define classical behavior. 
Equations of motion predict correlations in time, which in quantum mechanics are 
properties of the probabilities for time histories. The statement that the earth moves 
on a classical orbit is most correctly understood in quantum theory as the assertion 
that, among a decoherent set of coarse-grained histories of the earth's position in time, 
the probability is high for histories exhibiting the deterministic correlations in time 
implied by Newton's laws and low for all others. To discuss classical predictability 
therefore we should be dealing with the probabilities of time histories, not merely with 
the time dependence of expected values. 
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• The Ehrenfest derivation relies on a close connection between the equations of motion 
arising from the fundamental action and the phenomenological equations of motion 
determining classical correlations. There is no such connection in general. In general 
situations we expect classical equations of motion like the Navier-Stokes equation re- 
lating values of continuum hydrodynamic variables at different times, incorporating 
phenomenological equations of state, and exhibiting dissipation, noise, and irreversibil- 
ity. The equations of the fundamental theory, whether one takes it to be quantum 
electrodynamics or M-theory, exhibit none of these phenomena and are at best only 
distantly related in form. We need to derive the form of the classical equations of 
motion as well as the probabilities with which they are satisfied. 

• The Ehrenfest derivation posits the variable — the position x — in which classical de- 
terministic behavior is exhibited. But the quantum mechanics of any closed system 
will exhibit a great many complementary sets of decoherent histories some of which 
may exhibit deterministic correlations in time. What distinguishes coarse-grainings in 
terms of the familiar quasiclassical variables from all other possibilities exhibited by 
a closed quantum mechanical system? Certainly it is not their relation to the vari- 
ables of the fundamental theory, which is typically only distant, as described above. 
Rather, it must lie in the relative utility of different coarse-grainings for prediction, 
with quasiclassical variables being highly predictable. A complete derivation of classi- 
cal predictability must seek to distinguish classical coarse-grainings from all others. 

• The Ehrenfest derivation deals with the expected outcomes of "measurements" on 
an otherwise isolated subsystem. However, in quantum cosmology we are interested 
in classical behavior in much more general situations, over cosmological stretches of 
space and time and over a wide range of systems including the universe as a whole, 
whether or not they are receiving the attention of observers. We are interested in sets 
of alternative histories that can be assigned probabilities whether or not they describe 
measurement situations; in the quantum mechanics of closed systems that means the 
sets must decohere. Decoherence is thus a prerequisite for classical behavior. 

Histories of the quasiclassical domain of everyday experience are coarse-grained by values 
of usual quasiclassical variables. These include various sorts of hydrodynamic variables — 
averages over suitably small volumes of densities of conserved or nearly conserved quantities. 
Densities of energy, momentum, baryon number, nuclear and chemical species are examples. 
The system behaves classically when the probability is high for histories that exhibit corre- 
lations in time summarized by phenomenological classical equations of motion such as the 
Navier-Stokes equation. 

Simple arguments suggest why histories of these hydrodynamic, quasiclassical vari- 
ables should decohere and exhibit classical correlations in time. Coarse- graining is generally 
necessary for decoherence. Roughly speaking, a coarse-graining divides the variables of the 
system into those that are followed by the coarse-graining and those that are ignored. The 
ignored variables constitute the environment. Interaction between these classes of variables 
is necessary to dissipate the phases between different coarse-grained histories and to achieve 
decoherence. However, that same interaction produces noise, which causes deviations from 
predictability. Integrals of densities of conserved or approximately conserved quantities are 
natural candidates for quasiclassical variables. Their approximate conservation enables them 
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to resist deviations from predictability caused by the noise arising from their interactions 
with the rest of the universe. Further, following standard arguments of nonequilibrium sta- 
tistical mechanics, their approximate conservation leads to correlations in time summarized 
by a closed set of equations of motion. All isolated systems approach equilibrium. However, 
averages of approximately conserved quantities over suitable volumes approach equilibrium 
slowly. Closed sets of equations of motion result when the volumes can be chosen large 
enough that statistical fluctuations and noise are unimportant, but small enough that equi- 
librium is established within each volume in a time short compared to the dynamical time 
scales on which the variables vary (see, e.g. ||). The constitutive relations defining equa- 
tions of state, coefficients of viscosity, diffusion, etc., are then defined, permitting closure of 
a set of hydrodynamic equations of motion. Local equilibrium being thus established, the 
further equilibration of the volumes among themselves is governed by these equations. 

Despite the plausibility of the above general and simple qualitative picture, its validity 
has been only partially investigated quantitatively. To make the argument quantitative, in 
light of our earlier discussion, requires at a minimum an investigation of sets of histories 
coarse-grained by ranges of quasiclassical variables which has the following features: 

• Establishes the decoherence of sets of histories sufficiently coarse-grained by ranges of 
quasiclassical variables. 

• Establishes with high probability deterministic correlations in time summarized by 
closed systems of classical equations of motion for reasonably realistic initial conditions. 

• Compares the decoherence and predictability of different coarse-grained sets of histo- 
ries, both within the family of coarse-grainings by quasiclassical variables and with 
other coarse-grainings of distinct character. 



Despite intensive investigation of all of these points separately there is as yet no analysis 
which combines all three. There are many investigations of the mechanisms of decoherence 
of histories |IO|,^j5|Jll[] and of the closely related decoherence of density matrices |12[ , espe- 



cially in linear systems. However, these studies have typically posited a fixed division of the 
fundamental variables into those describing a "system" and those describing its "environ- 
ment" . Coarse-grainings follow variables of the system while ignoring those constituting the 
environment. Such a fixed system-environment split is intuitively accessible and correctly 
models many mechanisms of decoherence, but it is not general. 

Rather, coarse-graining is the general notion which, when possible, determines a family 
of system-environment splits. Different coarse-grainings lead to different possible notions of 
system and environment, that division is not usually unique, and for some kinds of coarse- 
graining no system-environment split is possible at all. (See Appendix A.) Even when a 
system-environment split is possible at one time, different system-environment splits could 
be needed at other moments of time. A fixed system-environment split is therefore neither 
general nor necessary for formulation quantum mechanics. However, as such workers as 
Feynman and Vernon |I3| |, Joos and Zeh |14j], Zurek ||15|| , Caldeira and Leggett ||16|| , and 
Omnes |17| fully appreciated, a system-environment split is an important tool for analyzing 
specific coarse-grainings and for understanding the physical mechanisms of decoherence. We 
shall utilize this tool extensively is what follows. 
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The emergence of deterministic correlations in time governed by classical equations of 
motion has been investigated for fixed system-environment splits @|. In a recent elegant 
paper, Halliwell p] has derived classical equations of motion for hydrodynamic variables 
although with a special assumption about the nature of the environment. Neither of these 
works compared quasiclassical coarse-grainings with nonclassical ones. In their pioneering 
paper on classical behavior in systems of interacting spins, Brun and Halliwell compared 
a family of coarse-grainings but did not derive classical equations of motion. 

This paper moves the analysis of decoherence and classicality a step towards realis- 
tic coarse-grainings by hydrodynamic variables in the context of a simple model — a one- 
dimensional chain of point masses with linear, nearest-neighbor interactions. (See Figure 
1). This is the simplest model which exhibits time correlations governed by a continuum 
equation of motion — the wave equation — in appropriate coarse-grainings for certain initial 
states. Yet the Lagrangian of the model is quadratic in all coordinates, so it can be tractably 
analyzed with standard path integral techniques [T^JT^J^]. Although our analysis is only for 
this very simple system, and does not deal with coarse-grainings by hydrodynamic variables 
per se, it does display all three features listed above. 

The details of the model are laid out in Section II, where we introduce a family of coarse- 
grainings. In the simplest, the chain of M identical particles is divided into Ai groups, 
each consisting of iV neighboring atoms. Histories are partitioned by ranges of values of 
the center of mass displacement in each group. We find that the effective equations of 
motion for these variables are well approximated by the classical wave equation. We obtain 
quantitative estimates for the statistical noise causing deviations from the predictions of the 
wave equations. 

We compare this quasiclassical coarse-graining with other members of a family of sets of 
alternative coarse-grained histories constructed as shown in Figure 2. The chain is divided 
up into M. groups of N atoms, each group consisting of equally spaced clumps of d atoms. 
We coarse-grain by equal ranges A of the values of the average (center of mass) positions of 
the atoms in a group at equally spaced intervals of time At. A family of sets of alternative 
coarse-grained histories is thereby defined parameterized by (N, d, A, At). When d = N the 
N atoms in a group are all neighbors, the coarse-graining entirely local, and the average 
position related to the approximately conserved center of mass momentum of the group. 
There is nothing special about the remaining members of the family except that they are 
amenable to analysis and range from local to highly nonlocal as d decreases from iV to 1. 

We find that for each member of this family that decoheres the probabilities of the 
histories can be thought of as obeying classical equations of motion augmented by noise [cf 
(O)]. We compare the members of this family with respect to three properties bearing on 



classical behavior: decoherence, noise, and the computational complexity of their equations 
of motion. We find marked differences which are describe in detail in the Conclusion, 
but which we summarize briefly here. As measured by the smallness of the ratio of the 
decoherence to dynamical time scales, decoherence is not a major constraint on predictability 
for "realistic" coarse-grainings. Noise interferes with the predictability of the nonlocal coarse- 
grainings much more than for the local ones. The equations of motion of the local coarse- 
grainings require many fewer steps to evolve over a given time to a given accuracy than 
do those for the nonlocal coarse-grainings. In short, we find that the familiar local coarse- 
grainings are more predictable than the nonlocal ones. 
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II. THE LINEAR CHAIN 

A. The Chain and Its Coarse-Grainings 

In this section we lay out the details of the model. We consider a linear chain of M 
"atoms", each of mass /i, separated at rest by equal distances Ax, and with displacements 
Xi, i = 0, ■ • • , M — 1. (See Figure 1.) Each atom interacts linearly with its nearest neighbors. 
The action describing this system is: 

S[x j (t)}= 1 -^Y: f dt {(x,(t)) 2 - ^ [x j+1 (t) - x 3 (t)] 2 } (2.1) 

/ j=0 J ° 

where, using periodic boundary conditions, we take x_\f = Xq. 

The linear harmonic chain is the simplest model of a solid and is treated in many refer- 
ences.]] The Xj(t) may be conveniently analyzed in terms of spatial modes 

Af/2 

x j( t ) = E [at(t)ft(j)+c.c] 

1=0 

where 

ftU) =^ 2 exp(2mj£/Ar) . (2.3) 
These are the normal modes of the chain and the ae(t) oscillate with frequencies 

u)i = 2ui sm(7i£/J\f) 
The action expressed in terms of normal modes takes the simple form 

S[at(T)) = -v'£ dt (\a t {t)\ 2 ~ ^W{t)\ 2 ) ■ (2-5) 

Both classically and quantum mechanically we will compare members of a family of 
coarse-grainings defined in part by two integer parameters N and d. The M particles in 
the chain are divided into A4 — J\f/N groups of N atoms each. Each group consists of N/d 
clumps of d particles spaced by M.d. (See Figure 2.) When d = N there is only one clump. 
This coarse-graining is local — the atoms are as closely spaced as possible in the chain. When 
d — 1 there is only one particle per clump. This coarse-graining is nonlocal, as the atoms 
are as spread out as possible in the chain. The variables followed by the coarse-graining are 
the average displacements of the groups. 

^ N/d-l [d/2] 

X f(t)=K E E XJd+m+kMd(t) ■ (2.6) 

k=0 m=[-d/2]+l 



(2.2) 



(2.4) 



1 For a classic reference see, [18]. 
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FIG. 1. The microscopic picture of the one-dimensional harmonic chain. N equal masses (i are 
spaced an equilibrium distance Ax apart. The displacement of the jth mass is denoted Xj. The 
masses have a linear restoring force between them with a spring constant k = [jluj 2 . The chain is 
assumed to have periodic boundary conditions. 



The first sum is over the clumps in the group, and the second sum is over the atoms in 
a clump. The range of the sum over m depends on whether d is even or odd. For odd d 
it ranges from —{d — l)/2 to {d — l)/2. For even d the range is —d/2 + 1 to d/2. The 
brackets [x] in ( |2.6|) denote rounding down to the next lower integer. The range of J is 
J — 0, • • • , M. — 1 making a total of M. coarse-grained coordinates. 

We consider sets of histories constructed from exhaustive sets of exclusive equal ranges 
A of Xj® separated in by time intervals At. There is therefore a family of coarse-grainings 
each member of which is labeled by values of the four parameters (N, d, A, At). 

It is convenient to expand the coarse-grained coordinates in terms of spatial modes as 
well 

M/2 

Xf\t) = £ [Af{t)F^\j)+c.c] (2.7) 

L=0 

where 

F { L d) (J) = M- 1/2 exp(2niJL/M) . (2.8) 

The parameter L labels the coarse-grained mode just as i labels the fine-grained mode 
at. The corresponding frequency depends on the choice of coarse-graining 
will be a superposition of the normal modes of the chain. Specifically, 

Af{t) = Y^cf k a m{k)N/d {t) (2.9) 

k=0 

where 

, M (L + kM/2, k even , , . 

mik) = \-L+(k + l)M/2, k odd . (2 ' 10) 

The coefficients cik are given by 
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N=6 



d=6cccxxn 



d=3 ooa 



d=2 oa 



d=i 



FIG. 2. The family of coarse-grainings under consideration. The total number of masses Af 
is divided into Ai groups of N masses each. These groups are then further subdivided into N/d 
clumps of d masses each. The clumps of a given group are spaced so that clumps of all the other 
groups occur before a clump is repeated, that is, by a distance Md. The coarse-grained variable 
Xj is the average displacement of the masses in the Jth group. In the figure, N = 30, Ai = 5, 
and N = Q. Masses in the same group are labeled with the same symbol, with a different symbol 
(including the shading) for each group. The arrangements are shown for d = 1,2, 3, 6, with d = 1 
being completely delocalized and d = 6 being completely localized. 



1 sin [irm(k)/M] 



J im(k)Lp{d) 



d^N sm[irm(k)/(Md)} 



[2.11) 



with the phase ip(d) vanishing if d is odd and equal to — n/ (Aid) if d is even. 
In the delocalized limit, d—1, the are just single normal modes 

A { l\t) <xa LN (t). 



(2.12) 



By contrast, in the localized limit d = N, the A^\t) are superpositions of all modes higher 
than L, viz. 



N-l 

A[ N \t)= £ cf k a m{k) (t). 



(2.13) 



fc=0 



These facts will be useful in what follows. 



B. Environments 



For several purposes in the subsequent analysis it will be convenient to divide the config- 
uration space of fine-grained modes {a^} up into a subspace spanned by variables that are 
followed by the coarse-grainings introduced above, and and a subspace spanned by variables 
that are ignored. The M. variables {A^} may be said to define the configuration space of the 
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"system". There is considerable latitude in the choice of variables {qjfi}, a = 1,2, ■ ■ ■ ,N—M. 
that define the "environment". The only requirement is that the {^4^} and the {q^} span 
the whole configuration space of the fine-grained modes {ag}. A convenient choice for our 
calculations will be to take the to be all but the lowest fine-grained modes that con- 
tribute to the through That is, we take the qf$ to be the real and imaginary 
parts of: 

a m (k)N/d, k=l,2,---,d-l; L = 0,- ■ ■ ,M/2 . (2.14) 

Each label a corresponds to a pair (L,k); the exact form of the correspondence will not 
be important for us. Thus defined, the environment variables are normal modes which 
classically obey: 

q a d) + U d) ) 2 q a d) =0. (2.15) 

The definition Q2.14| ) defines a different environment for each d, and even for a given d it is 
but one of many possible choices, though the effect on the system must be independent of 
this choice. 

Since only d fine-grained modes contribute to each coarse-grained mode A^\ this "system 
plus environment" only includes the full set of modes in the localized case d = N. For all 
other d, there will be some fine-grained modes that contribute to none of the coarse-grained 
modes, and which thus can be neglected. When d = 1, each coarse-grained mode 
depends on only a single fine-grained mode, and hence has no environment at all. 



C. The Initial State 

As we discussed in the Introduction, the classical behavior of any quantum system is 
exhibited through the probabilities of certain decoherent sets of alternative coarse-grained 
histories arising from particular initial conditions. In this paper we compare the predictabil- 
ity of the different members of the family of coarse-grainings introduced above for a single 
class of initial conditions. That is like the situation in cosmology, where the initial condition 
is given by fundamental theory and we search among the possible sets of coarse-grained 
alternative histories for sets that exhibit more predictability than other sets. 

In this simple model of the harmonic chain there are no interactions between different 
normal modes. This is an idealization; realistically there will be interactions. These inter- 
actions mean that the normal modes will tend to equilibrate. The highest frequency modes 
can be expected to come into equilibrium first, the lowest frequency modes last. Thus, at 
an intermediate time in the approach to equilibrium, the chain can be described by a state 
in which a range of low frequency normal modes are not yet in equilibrium, but with higher 
frequency modes that are in equilibrium at a temperature T. 

Indeed, a garden variety string would be described in this way. In a plucked violin string, 
for instance, the wavelengths that are reasonable fractions of the length of the string are not 
in equilibrium, but very short wavelengths — short compared to the transverse dimensions of 
the string, certainly — are likely to be in equilibrium. 

Classically an initial condition is given by a distribution on the phase-space of mode 
amplitudes ag and their conjugate momenta Tig. We denote these collectively by {zg} = 
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{(a£, He)}, t = 1, 2, • • -TV. We assume that the significant expected values of Z£, denoted zp, 
occur only low frequencies, so that for I greater than a cutoff £c we have z# ~ 0. We assume 
thermal fluctuations around these expected value. A distribution which gives this is: 

P(M) = nWe-*"* - *^) • (2-16) 

Here Hp(zp) is the Hamiltonian (n]/fi + fiujaj)/2 for each mode, and Zp is a normalizing 
factor. 

The corresponding assumption in quantum mechanics would be a density matrix of the 
form: 

JV-l 

p= H exp[i(d* e 7c e - TT* e a e + h.c.)/2] {Zj x e~ Hl l T ) exp[-i(a}ft t - n* e a e + h.c.)/2] (2.17) 

1=0 

where hp and Tip are the quantum operators corresponding to the classical variables zp = 
(ap, Tip), and the exponentials sandwiching the thermal state are phase-space displacement 
operators. These are the initial conditions we assume in our analysis. 

A system is said to be in thermal equilibrium when its state is such that the probabilities 
of quasiclassical alternatives are accurately reproduced by a thermal density matrix. Typical 
examples are coarse-grainings by volume averages of densities of conserved quantities such 
as energy, momentum, etc. However the probabilities of an arbitrary set of alternatives will 
not be generally be reproduced by the thermal density matrix. It is plausible that many 
(or perhaps most) initial states of the linear chain relax after a suitable time to a state of 
local thermal equilibrium for which the probabilities of the local coarse-grainings described 
above are reproduced by a density matrix of the form ( |2.16| ) or ( |2.17| ). However, we have 



assumed something much stronger which is that the probabilities of the whole family of 
coarse grainings are accurately reproduced by these density matrices. That is a much more 
restrictive condition on the initial condition. 



III. EQUATIONS OF MOTION FOR CLASSICAL COARSE-GRAININGS 

In this section we consider the classical linear chain. Each atom moves according to 
a classical equation of motion. However, with a probabilistic initial condition, a set of 
coarse-grained variables (such as the Xj discussed in the previous section) need not obey 
any closed system of deterministic equations. We analyze the constraints on the initial 
distribution p(z°) necessary for the {Xj} for each d to obey and closed set of deterministic 
equations and we exhibit the form these equations take. The general kind of analysis we 
describe has been considered by many authors, for example Zwanzig |H| and Brun ||. What 
is new here is the application of these methods to comparing different coarse-grainings of 
the linear chain. 



A. The Probability of Determinism 

A set of functionals of the phase-space paths FA[z(t)], A = 1, • ■ ■ , P and a set of ranges 
{A a } in R p define the most general kind of classical coarse-graining. The probability p a 
that the functionals Fa have values in the range A a is 
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Pa = / 6ze a (F A [z(t)]) 5[z(t)-z t (z°)}p(z ) (3.1) 



Here, e a is the characteristic function (1 inside, outside) for the range A a and z t (z°) is the 
classical evolution of the initial data z°. The functional 5-function enforces that evolution, 
assigning zero probability to all paths which do not conform to it. The functional integral 
is over all phase-space paths z(t) including an integral over the initial conditions z°. 

Utilizing this framework one can calculate the probability that the coarse-grained vari- 
ables such as the {Xj } exhibit deterministic correlations in time. The functionals Fa for 
example might defined different orbits in phase-space. Even more simply, one can evaluate 
the probability that a set of deterministic equations of motion^ 

Sj{t,X(r)]=0 (3.2) 

holds for a set of variables {Xj®} (denoted without indices in ( |3.2|) ) at time t, by calculating 
the probability that the functionals Sj(t, X(t)] have the value zero. If that probability is 
high then the coarse-grained variables obey the deterministic set of equations 



B. Equations of Motion 

We now consider equations of motion for the coarse-grained average positions xf\t) 
introduced in Section II. Equivalently, and more conveniently, we can consider the equations 
of motion for the corresponding spatial modes (t) defined by ( |2.7| ) and given terms of the 
fine-grained modes by (|2.9|) . It is convenient to take the system-environment split defined 
by ( p.!4| ) in which case the ignored coordinates are the {q^} of ( |2.15| ). 



The normal modes of the chain ag oscillate with the frequencies ue given in eq Q2.4|) . 
Evidently as a consequence of (|2[ 



d-l 



A L{t) — - °Lk ^m(k)N/d am(k)N/d(t) . (3.3) 



fc=0 



We have suppressed the superscript (d) in this equation, as we shall do for clarity in the 
remainder of this section. The ignored coordinates were chosen to coincide with all but the 
lowest frequency normal mode in (|3.3|). Therefore that equation may be rewritten as 



A L (t) + n\A L {t) = f L {t) , (3.4) 

where 

= ULN/d , (3.5) 

and 



2 We follow the convention that f(t, x{t)\ means that / is a function of t and a functional of X(r) 
for example Jq x(r)dr. 
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h(t) = J2(nl-ul)c La q a (t) . (3.6) 



a=l 



Here, u a , ci a , etc. are the values of ue, clu appropriate to the value a. 

The simple equation of motion (|2.15 ) can be solved to express the q a {t) in terms of their 



initial values 

p a (0) 

q a {t) = g o (0) cos (u a t) + — — sin (u a t) . (3.7) 

The q a (t) are probabilistically distributed with a distribution that depends on the distribu- 
tion of the initial conditions through (|3.7|) . Let Fi(t) be the expected value of /z,(t), 



F L (t) = M[f L (t)} (3.8) 

where M[/^] denotes the mean over the distribution of initial conditions. (We use M[-] 
for classical expected values and (■) for quantum mechanical ones.) Then equation Q3.4| ) 
becomes 



where 



S L (t, A l {t)\ = A L (t) + n\A L {t) - F L (t) = Af L (t) , (3.9) 



Af L (t) = f L (t)-M[f L (t)} . (3.10) 



We can think of S^t, A L (r)] = as a deterministic equation of motion and Af L as 
classical stochastic noise causing deviations from that determinism. Since the expected 
value of £l(1;,Al(t)] is zero, the coarse-graining will approximately obey the deterministic 
equation £l(£, Al(t)] = provided the noise is small. A simple measure of its size is 

M[Af L (t)Af L (t')} . (3.11) 

The situation is even simpler if, in addition, i*x(t) is negligible, as it is in many realistic 
situations. Then the chain obeys the equation 

A L (t) + n 2 L A L (t) = Af L (t)^0 (3.12) 

with high probability — an equation whose form is independent of the initial distribution. 

For our chain model F^t) does not vanish. Because the initial condition is diagonal in 
the fine-grained modes rather than the A^s and q a y s, there is a non- vanishing contribution 
to the driving force Fl proportional to (Al(0) — ML4l(0)]). However, even in terms of At 
and q a , the initial condition is close to diagonal for large N and Ai; thus, this term is small 
compared both to the noise and the terms of the homogeneous equation of motion ( 3.12 ). 
The time average of Fi(t) also vanishes. Hence, we can safely neglect it. 

Unlike effective classical equations for many systems, eq ( |3.12| ) exhibits neither nonlo- 
cality or dissipation. These simplifying features can be traced to the conservation of the 
energies of the individual normal modes. The environment — the g a 's — are all normal modes. 
Energy, therefore, cannot be exchanged between the A L and the q a on average, that is be- 
tween the system and its environment. Defining the energy for a coarse-grained mode is 
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somewhat ambiguous, but simple definitions give an energy which fluctuates about a fixed 
average. This is special to this linear system which has a high level of integrability. 

The initial distribution p{zf) determines whether the classical noise is small. The noise 
will be negligible when there is a negligible probability for any significant initial excitation 
for the q a . Thus, for initial distributions that favor low frequency modes below some cut-off 
value, we expect low noise and deterministic behavior of the Ar,(i). We shall return to this 
question in more detail in Section V. 



C. Deriving the equations of motion from the Lagrangian 



While the above derivation is perfectly correct, it is in some sense a short-cut to a more 
standard procedure. In this, we would substitute the change of variables fl2.7| ) and fl2.14| ) 
into the expression for the action ( |2.1|) , and derive the Euler-Lagrange equations for the 



new variables. This is important in making contact later with the quantum case, since that 
derivation proceeds from the classical Lagrangian. 



The action (|2.1| ) may be expressed more compactly in matrix notation writing a(t) for 
the vector with components a^(i). Then 



S [a(r)] 



1 
2 



dt 



(a fid) — (d n 2 uj 2 d 



(3.13) 



We use a Hermitian inner product indicated only implicitly, so the (ab) = X^a|&g. Sand- 
wiched in inner products like those of ( [3.13D , uj 2 and \i denote the diagonal matrices of normal 
mode frequencies and masses (all equal) respectively. For this case, fj, 1 ^ 2 ^ 2 /x 1 / 2 



The coarse-grainings are by the variables Xj , or equivalently by their modes A£ J . The 
coordinates qjf* are the fine-grained modes that are not followed by the coarse-graining. 
Continuing to suppress superscript (d), we write 



i(d) 



at 



S eL A L + T eaQa 



(3.14) 



or in matrix notation 



SA + Tq 



(3.15) 



In general, there is considerable arbitrariness in the matrices S and T corresponding to 
how the variables q a are chosen. The only requirements are that ( |3.15| ) must reproduce the 
definition of the A L 's ( |2.6| ) and (|2.7]) and the set of (A L ,q a ) must span the space of the a£. 
The choice (|2.14j ) specifies a definite T and S. 

Inserting the transformation ( |3.15| ) into the action ( |3.13| ) one has 



S 



A(t),q(t) 



dt 



A M ss A) + [qM TT q 
+ (lM 5T <fJ + ($M TS A) - (AVssA 
- (qV TT q) - [AV ST qJ - (qV TS A 



(3.16) 
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where 

M AB = A f n B , (3.17a) 

and 

V AB = AV^V^B . (3.17b) 

With the action in this form, it is a simple matter to derive the equations of motion for 
A and q: 

M SS A + Mst'I = -VssA - V ST q , (3.18a) 

M TT <f + M TS A = -V TT q- V TS A . (3.18b) 

Eq. ( |3.18| ) looks complicated, but in fact we can easily recover our earlier result. Because 
the q variables are normal modes, they obey the usual harmonic oscillator motion ( 2.15 ). 
By plugging the second equation into the first, the equation becomes 

(M ss - MstM£M ts )A+ (V ss - M ST M£V TS )A = -(V ST - M 5T M^V TT )g. (3.19) 



If we then plug in the simple harmonic solution ( |3.7|) for the ignored modes q the equation 
reduces to the simple Q3.12p above multiplied by the diagonal matrix (M S5 — Ms T M T ^M TS '). 
(See Appendix B for the proof of this.) 

This is a slightly more complicated derivation of the equation of motion, but it still relies 
on the shortcut of knowing the solution for the motion of q. That will not be available in 
quantum mechanics. Without making use of this, we could solve the linear equation (|3.18b[) 
above for q(t) and insert it into (|3.19| ). The solution for q(t) is 

M^qXt) = cos[f2t]M^g(0) + fT 1 sin[Ot]M^(f(0) 

-ft 1 C dt' sin[Q(t - t')]M^ /2 {M TS A{t') + V TS A(t')), (3.20) 
Jo 

where 

O 2 = M^ /2 V TT M^ /2 , (3.21) 
and it gives us the very complicated-looking equation of motion 

(M 3S - MstM£M T s)A+ (V 5 s - M ST M^V TS )A 

/•t 

+ {Vst - MstM^Vt^M^O" 1 / dt' sin[ft(t - t')]M r y /2 (M TS i(t') + V TS A(t')), 

Jo 

= -(V ST - M 5T M^V Tr )M^ /2 (cos[Ot]M^g(0) + IT 1 sin[Ot]M^^(0)). (3.22) 

Two observations can be made. First, the solution ( p.20|) while still harmonic motion, 
appears to be at quite different frequencies from the form ( |3.7| ), and also includes a driving 
term absent in that case. Second, equation ( |3.22|) appears very different from the rather 
simple harmonic oscillator equation ( |3 . 1 2| ) . These differences, however, must be apparent 
rather than real; both ( |3.7|) and ( 3.20|) follow from the same Lagrangian, as do ( P-12 ) and 
( ^.22| ). In fact, the more complicated form ( |3.22| ) of the equation of motion is related to the 
simple form ( |3.12| ) by an invertible transformation, and thus has exactly the same solutions. 
(See Appendix C.) 
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D. Transforming the classical noise 



The fact that classical equations of motion can be represented in widely different forms 
is nothing new, of course. But the apparent complexity masking the comparatively simple 
underlying dynamics complicates the analysis of the quantum case. For the classical case, 
it is best to stick with the simple form fl3.12| ). 



The two forms do differ in one important respect. The retarded term in equation ( 3.22 



is most readily identified as part of the homogeneous equation of motion, while in ( J3.12 ) it 



is implicitly included in the noise terms. We shall see that this ambiguity in the definition of 
the noise becomes important in comparing the classical and quantum cases: one definition 
is most natural in the classical case, while the other is most natural in the quantum case 
(where A(t) need not obey the classical equation of motion). 

Let us restrict ourselves to a single coarse-grained mode A^f (t) and its associated modes 
q(t). If we define the two forms of the noise to be the right hand sides of equations ( |3.12| ) 
and (|3.22| ), respectively, they can be expressed 

Af L (t) = c - U 2 Q ) (cos[n Q t]Aq(0) + tl^ sin[0 Q t] A'q(0)) , (3.23a) 

and 

Af' L (t) = Nc(n 2 L I - H|) M^ /2 (cos[Ot]M^A(f(0) + fT 1 sin[Ot]M^A^(0)), (3.23b) 

where c is a vector in the space of the g's with elements equal to the coefficients c&, Oq is 
a diagonal matrix on the space of the g's with diagonal elements Ub, and ft is the effective 
frequency matrix defined by ( |3.21|) . The structure of these two expressions is closely parallel, 



and in Appendix C we show that we can switch from one form of the noise to the other by 
means of an invertible linear transformation. As far as determining the classical dynamics 
and predictability, they are equivalent. 

However, if we look at the absolute strength of the noise as a function of L or d, the form 
of the noise can make a difference. We assess this by looking at the correlation function 
( [3. 11| ) for the two forms of the noise. 

In the initial state (|2.16j ) we find the expectation values 

M[Ag b (0)Ag*,(0)] = k B T{V£) w , 

M[Aq b (0)Aq*(0)] = k B T(M£) w . (3.24) 
We can use these to calculate the correlation functions ( [3.1 H ) 

M [Af L (t)Af L (t')] = k B Tc(n 2 L l - H|) (cosfOQtjV^cosfOc/] 

+fi Q 1 sm[n Q t]M^ sinfOgt^g 1 \ - H|) c , (3.25) 

and 

M [Af' L (t)Af' L (t')} = Nk B Tc(p? L l - n 2 Q ) M^ 2 (cos[nt}n- 2 cos[Ot'] 

+fT x sin[Ot] sin^O'MM^ 72 (Q 2 L I - c . (3.26) 
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FIG. 3. The time-averaged mean force squared of the noise, S 2 , for the harmonic chain as 
a function of the coarse-graining d. This noise is chosen for three typical coarse-grained modes 
L = 30, 65, 100, with M. = 630 and TV" assumed to be very large; N dependence is absorbed into 
the units, /cgTa; 2 /Nfi. Note that the noise vanishes for d = 1, but otherwise assumes its highest 
values at low d, dropping off like 1/d at high d, as is predicted by the analytical result ( |5.6| ) . 



Using these expressions, we can estimate the mean strength of the noise by taking the 
average over time, 

S 2 = lim I f 1 M[(f L (t) - F L {t)f]dt . (3.27) 

In Fig. 3 we plot S 2 as a function of the coarse-graining d. We see that the noise strength 
falls off steeply as a function of d. Thus, the noise becomes lower in the highly localized 
case, and the motion of the localized coarse-graining is most predictable. Later we shall see 
that the Lagrangian form of the noise is closely related to the strength of decoherence in the 
quantum case, and closely resembles this "simple" noise. 
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IV. THE WAVE EQUATION 



The classical wave equation for the string does not follow directly from the deterministic 



equations for the chain (|3.12|) . A restriction on the initial distribution is required beyond that 
necessary for determinism. This is the requirement that only very long wavelength modes 
of the chain are excited. We assumed such a restriction on the initial distribution in Section 
IIC, but in this section we will examine the requirements on the short wavelength cutoff £q- 
This small gradient approximation is necessary for the validity of the wave equation as for 
other familiar continuum equations such as the Navier-Stokes equation.^ 

A. The Small Gradient Approximation 

The expansion (|2.7| ) of the coarse-grained coordinates Xj® (t) in terms of coarse-grained 
modes A^\t) connects the equation of motion ( |3.12|) for the A^\t) to an equation of 
motion for the Xj d \t). The character of these equations for the Xj d \t) is determined by 
the dispersion relation for which is [cf. ( p.5|) , ( 2.4) ]: 



4 d) =2-sin(-^) . (4.1) 



Since ( |4.1|) is not linear in L, the equations for xf\t) following from ( [3 . 1 2|) will generally 
relate XfP(t) to all the other coordinates Xj^(t) of the chain. The form of the equations 
for Xj d \t) simplifies if the initial distribution is such that only modes with 

£<Af (4.2) 

have any significant probabilities, that is £ c <C N '. Then from fl2.4|) , uji ~ 2iruj£/J\f and in 
particular 

{d) 2nuL 

Md • {A - 6} 

In this small gradient approximation, the equation of motion for the Xj d \t) implied by 
( CT ) and (g/7|) is 



Xf{t) = ^ [4t(t) - 2Xf\t) + Xf_ x {t)] ■ (4-4) 

Only nearest neighbor interactions are involved in the small gradient approximations. 

Usually a much stronger condition is meant by the small gradient approximation, namely 
that the only modes with significant probabilities are those with 

£<A^<AT, (4.5) 



5 See, e.g. (20). 
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that is £c <C M.- For d = N this condition ensures that many groups will fit into a 
wavelength so that the Xj vary only slightly from one J to the next. We will see below that 
this is essential to deriving the continuum wave equations. 

Condition ( f4.5|) can be very much stronger than (|4.2|) . In a 10 cm length of string with 
typical interatomic spacings M ~ 10 9 . Dividing the string into .1 mm lengths constituting 
the groups in the d = N case gives N ~ 10 6 . 

The condition ( |4.5| ) would imply that no values of [t] would be excited above L = 
since the i values contributing to any d = 1 mode are all larger than N [cf. (|2.9| )1 except 
i = 0. That makes the equation of motion ( |4.4| ) in the delocalized case not incorrect, but 
rather trivial since both sides are negligible. 



B. The Continuum Approximation 



For the localized coarse-graining d = N, the difference equation ( |4.4| ) is well approxi- 
mated by the wave equation when N is large. The derivation is standard, but we briefly 
repeat its essential features here. 

Recall that we denoted the mass of an atom by /i, and the spacing between atoms in the 
unexcited chain by Ax. The mass density a is therefore 



Multiplying both sides of 

aXf\t) = ^ 2 Ax 



a = fi/Ax . 
by a it can be written 

-x%\( t )-2xy\t)+xn(t) 



(NAxY 



(4.6) 



(4.7) 



The length NAx is the distance between the centers of the groups in the unexcited string. 
Assuming the strong form of the slow approximation gradient ( f4.5| ), we can approximate 
( [17D by the continuum equation 



d 2 X(x,t) ^d 2 X(x,t) 
c ^-^ = Y - 



dt 2 



dx 2 



(4.8) 



where Young's modulus Y is \ioj 2 Ax. Eq (|4.8|) is the wave equation for the propagation of 
compressional modes along the string. 



V. CLASSICAL PREDICTABILITY 

The fine-grained variables of a classical system obey a closed system of deterministic 
equations of motion. However, there is no guarantee that coarse-grained variables will. 
There may be no deterministic equations at all, or the set of those that do hold may not 
close, in the sense that there may not be a complete set of equations to solve for all the 
variables. In Section IIIB we demonstrated that the probabilities of the evolution of the 
coarse-grained modes Aj/(t) of the linear harmonic chain are reproduced by a classical 
equation of motion (|3.9| ) modified by noise: 
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Af{t) - & L Af(t) - Fi d \t) = Afi d \t) (5.1) 

for the family of coarse-grainings under consideration. The time evolution of the (t) will 
be classically predictable by the left hand side of ( |5.1|) if the noise term on the right hand 
side is negligible. In this section we analyze this requirement for classical predictability as 
a function of d. 

A simple estimate of the noise can be obtained by assuming a completely thermal initial 
state characterized by a temperature T. That is, we assume fl2.16| ) with Iq = 0. In this initial 
state, M[fi d) (t)} w 0. A measure of the magnitude of the noise fluctuations is M[(f {d) (t)) 2 ]. 
In a thermal state 

M[a t (t)a t ,(t)] = 5 u >k B T/(2iMJ 2 e ) . (5.2) 
We therefore have from (|3.6| ): 



M[(/f (t)) 2 ] = ^ £ \c Lk \^ {k) ( 1 - ^ ) f (u£ (Jk) V£ + M&) (5.3) 

v k=i V z 

where the final factor (arising from the initial condition) is very close to 1 for M. 3> L, that 
is, (M/2)(w| (fc) V^. + M^) kk = 1 + 0(L 2 /M 2 ), where £(k) = m(k)N/d is 

i(k) = (L + kM /2)N/d , for k even, (5.4a) 
£(k) = (-L + (k + l)M/2)N/d , for k odd . (5.4b) 

Every term on the right hand side of ( |5.3| ) depends on d although we have not indicated the 
dependence explicitly. 

Figure 3 shows a plot of the expected value of the square of the fluctuations (|5.3| ) as 
a function of d for some representative values of L. The plot is for A4 = 630; the value 
of iV cancels out of the plotted function, but is assumed to be very large in order to have 
sufficient factors d. As is evident from (|5.3| ), the noise vanishes for d — 1. That is because, 
as we noted in ( |2.12|) , is a single fine-grained mode not coupled to any other fine-grained 
modes. This special situation results from our idealized linear model of the chain. Even a 
tiny amount of nonlinearity would couple this mode to others and produce noise. 

Figure 3 shows that the noise generally increases with L and decreases with d. Analytic 
estimates can be derived from ( |5.3j ) when L ^ A4 < N. Assuming that the final factor in 
Q5.3D is well approximated by unity, we have for d = 2 

M[(/£° W) 2 ] - k -§f(nL/M) 2 , (5.5) 
while for large d comparable to 

M[(/r«« « ■ (5-6) 

The noise thus decreases inversely with d. For the "realistic" coarse-graining discussed at 
the end of Section IVA, where ~ 10 6 , the noise is 10 6 times smaller for d = N than it is 
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for d = 2. That is a vast advantage in predictability of the local coarse-grainings over the 
nonlocal ones. The origin of this advantage can be traced to the approximate conservation 
of the center of mass momentum of the local groups of atoms. 

The noise term on the right hand side of the equation of motion ( |5.1|) must be compared 
with the characteristic size of the deterministic terms on the left hand side to get a true 
estimate of the effect of noise on predictability. In ( 2.16 ) we assumed an initial state in which 
only fine-grained modes with I < lc had significant excitations above thermal noise. But, 
as we remarked earlier, the coarse-grained modes are superpositions of fine-grained 
modes with t > LN/d, as equation ( |2.9| ) shows. Thus, for d < LN/lc-, there will be no 
significant excitation of above the level of thermal noise. The subsequent dynamics 
is predictable only in a trivial sense. The string exhibits no motion except in response to 
thermal fluctuations. For true classical predictability d must be large, so the noise is low, 
but also so the deterministic terms in the equation of motion dominate the noise terms. We 
shall return to a more quantitative comparison of the noise and dynamical force terms in 
the Conclusion. 

Even leaving aside the competition between noise and predictability, the classical equa- 
tions of motion for nonlocal coarse-grainings are distinguished from the local ones by their 
computational complexity. Consider the classical equations ( |3.12|) for the coarse-grained 
modes, or equivalently, in the small gradient approximation, eqs. ( |4.4[ ) for the center of 
mass positions of the groups. The characteristic dynamical time scales are of order d/uj. 
The local coarse-grainings therefore vary the most slowly as would be expected from their 
association with approximately conserved quantities. Thus the computation of the evolution 
of the coarse-grained quantities A^ or Xj over a given time interval to a given accuracy 
will take a factor of N more time steps for the nonlocal equations than the local ones. For the 
"realistic" case when N ~ 10 6 that is a significant advantage in computational complexity 
for the local coarse-grainings. 

Thus, whether one considers the absolute value of the noise, the relative size of the noise 
and deterministic forces, or the effort needed to solve the classical equations of motion, the 
local coarse-grainings are more predictable than the nonlocal ones in the family we have 
considered. In the next section we turn to the same questions in quantum theory. 



VI. DECOHERENCE 

A. Quantum Mechanics of the Linear Chain 

In this section we turn to the classical behavior of the quantum linear harmonic chain. A 
quantum system behaves classically when the probability is high that coarse-grained histories 
exhibit the correlations implied by classical equations of motion. We will be concerned 
with the histories of the linear harmonic chain coarse-grained by values of the position 
averages {Xj(t)} defined in Section II. More specifically we shall consider, for each d, 
sets of histories defined by an exhaustive set of exclusive regions {A^ }, a fc = 1,2,3, ••• 

of the A^-dimensional configuration space spanned by the {Aj (t)} at a series of times 
tk, k — 1, 2, • • • , n with ti < t 2 ■ ■ ■ < t n . We take the same set of regions for each value 
of d and usually take them to be "cubes" in M. of equal sides A. An individual history is 
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then labeled by the particular sequence of intervals (a±, • • • , a&) = a. We then compare the 
probabilities that these histories are correlated by the classical equations of motion discussed 
in Section III for different values of d. 

Quantum interference between the individual members of a coarse grained set of histories 
must be negligible for probabilities to be consistently assigned to its individual members. 
When this condition is satisfied the set of histories is said to decohere. Decoherence is a 
prerequisite for probabilities. 

The interference between histories a and a' is measured by the decoherence functional 

D«V,a) = Tr • • • P d ai {t x ) pP^h) ■ ■ ■ P* n (t n )] . (6.1) 

Here, {Pa k {tk)} are an exhaustive set of mutually orthogonal Heisenberg picture projection 

operators projecting on the regions {A^ } of {Xj} at time The operator p is the 
Heisenberg picture density matrix of the system. The set of histories decoheres when the 
"off-diagonal" elements of are negligible, 

D {d \a',a) « , a' a . (6.2) 

There is an equivalent path integral expression for the decoherence functional of sets 
of histories coarse-grained by ranges of configuration space such as those under discussion 
here. Suppose the coarse-graining is entirely confined to times less than a final time tf, 
i.e., t n < tf. Then, for the linear chain coarse-grained by ranges of xf\t) 

D {d) (a',a) = J 5a! J 5a5(a' f -a f ) 

x exp (i {S [a{t')} - S [a(t)}} /%) p (4, oo) . (6.3) 

One integral in ( |6.3| ) is over paths a(r) on the interval [0,tf] which start at do at t = 
and end at a/ at t = tf including integrations over those end points. The integral is only 
over paths which pass through the regions {A^ fc } in {xj^} at the times The constraint 
on the Xj d ^ translates linearly into a constraint on the a through (|2.2j ). An integral over 



a'(r') is similar except that it is constrained by the coarse-grainings of a'. The sum in 
could have been expressed in terms of any configuration space variables. We have chosen 
the modes at because the action takes the simple form fl2.5|) . Equally well, we could have 
used the coordinates of the individual atoms, a^, % — 1, • • • , Af. 

Path integrals of the form ( |6.3|) have been extensively studied for quadratic actions and 
thermal density matrices by many authors [|l^,|n],|||| . All integrals are Gaussian in this 
situation and can be evaluated explicitly. The simplest way to review this is to recall a 
simple example. 



B. A Simple Example 

The sum-over-histories techniques used in this paper to calculate the decoherence func- 
tional for sets of alternative coarse-grained histories of the linear harmonic chain are gener- 
alizations of those used for more straightforward coarse-grainings of simpler linear systems 
In turn, these methods extend those pioneered by Feynman and Vernon and 
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Caldeira and Leggett flq| . While the results of this paper are algebraically more complex, 
the basic ideas are similar to those in simpler cases. To emphasize this connection with pre- 
vious work, and to explain the ideas in an algebraically simple context, we review a version 
of these simple models. 

The model consists of the linear chain under discussion but with a particle of mass M 
and position X(t) coupled by a linear interaction to atom 0. The total action is 

S tot [X(r),x{r)} = S[x(r)] + S x [X(r)] + S- mt [x (r),X(r)] , (6.4) 

where S [x(t)} is given by (|2.1j) , Sx[X(r)] is 

S x [X(t)] = f tf dr \ MX 2 (r) (6.5) 
Jo 2 

and the interaction term is 

S int [x q (t),X(t)] = -k I*' drX(T)x (r) . (6.6) 



We consider coarse grainings where only histories of X are followed. Thus the entire 
chain serves as the environment. 

The action 5[x(r)] is which is ( |2.5|) when expressed in terms of normal modes. 

The interaction is 

Sm [d(r),X{r] = -kAH r dra e (r)X(r) . (6.7) 



The problem summarized by ( |6.4j) can be mapped onto a problem studied by many au- 
thors following Caldeira and Leggett [I5|. We will follow the calculation of the decoherence 
functional described in j|] but in a notation that is designed to stress the analogy with the 
subjects of this paper. The translations are as follows, the first being the notation of this 
paper, the second of jO: 



a>e<-+ Qa , ui> <-> lo a , 

X(t) <r-> X(t) , [I *r* Wj , 

kM^ ^C a , M <-> M . (6.8) 

We assume an initial condition that is a product of a density matrix p(x' , Xq) for the 
particle and a thermal density matrix Pt{S' q , do) at temperature T for the environment. 

We consider a set of alternative coarse-grained histories defined by exhaustive sets of 
ranges for X, {A^} , {A^} • • • {A™} , a = 1, 2, • • • at a series of times t\ < ■ ■ • < t n . The 
decoherence functional for this coarse-grained set of histories is given generally by 

D(a', a) = J 5X' J 5X J 5a J 5a 5 (x' f - X f ) 5 (a' f - a)) 

x exp (z {S tot [X'(t) , a'(r)] - 5 tot [X{r) , a(r)]} /h) 

x P (Xq , X ) p T (4 , d ) . (6.9) 

The integral over the X's is restricted to paths that traverse the regions A^, , • • • , A" n defin- 
ing the coarse-grained history a = (a±, ■ ■ ■ ,a n ) and similarly for X"s. The integrals over 
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the a' e and ae are unrestricted and, given the assumed thermal initial condition, reduce to 
Gaussian integrals which can be carried out explicitly. The result, expressed in terms of the 
variables 

X(t)= l -[X\t)+X{t)] , (6.10a) 
£{t)=X'(t)-X{t) , (6.10b) 

£o * £o 



is 



D (a', a) = J a> 5x> J 5x e "t*M, W)~ p x + f , X - f (6.11) 



where 



A 



X{r) , £(r)] = -^M (dX/dt) | t=0 + /*' dt £(t)e(t , X(t)] 



+ ^ [ f dt[ f dt'at')ki(t'-t)m- (6-12) 
4 jo Jo 

The ingredients in eq. ( 6.12j) are the equation of motion 

e(t,X(r)] = -M^it) + J*dt'k R (t-t')X(t') , (6.13) 
together with the kernels 

k R {t) = -k 2 (fiAfy 1 ^^ 1 sin(u e t) , (6.14) 

i 

and 

kj{t) = ^(MO-^u^coth \^\ cos(^t) . (6.15) 

The imaginary term in ( |6.12j ) favors £(t) = 0, that is X'(t) = X(t). If it is large, the off- 
diagonal elements of D(a', a) will be negligible and decoherence achieved. The characteristic 

time scale tdocoh over which enough imaginary exponent is built up in between to make off- 
diagonal elements of D negligible is 

tdecoh ~ (kiA 2 ) 2 , (6.16) 



where ki and A denote characteristic sizes of the kernel ( |6.15| ) and the intervals A^,. 



If the set of histories decoheres, then the restrictions of the coarse graining on the inte- 
grals over £(t) in ( |6.11| ) defining the diagonal elements of D(a',a) can be ignored and the 



resulting Gaussian integral over £(£) carried out. The result is the following expression for 
the probabilities of the coarse-grained history a: 



p(a) = f 6X[det A; 7 /47rp 

J a 

x exp^~ dt' j\te(t',X(r)]kf v (t' - t) e(t, X(r))\ w(X , P ) . (6.17) 
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Here, kf v (t,t') is the inverse of kj(t,t r ) and w(X ,P ) is the Wigner distribution of the 
density matrix p. 

Eq. ( |6.17| ) shows that probabilities peak on histories obeying the deterministic equation 
of motion e(t,X(r)] = 0, but with a width in this distribution related to kj{t' — t). The 
same probabilities follow from a Langevin equation 

e(t,X(r)]+e(t) = , (6.18) 

with a stochastic noise force distributed with a correlation function 

M[£(t')£(t)] = hkj(t' - t) (6.19) 

Thus, decoherence and noise are connected. The stronger the coupling between system and 
environment the more rapidly interference between histories is dissipated, but also more 
noise which disturbs the deterministic dynamics and reduces predictability. 

The deterministic equation e(t,X(r)] = is exactly the same as the classical equation 
of motion which would be derived from the Lagrangian for this system, in the limit of zero 
noise. This result is generally true for systems with quadratic Lagrangians, and may hold 
approximately for systems with greater nonlinearities. More details than we have given here 
are presented in ||. We will follow the basic ideas of this model in our treatment of the 
family of coarse grainings of the linear chain. 



C. Decoherence Functional for the Linear Chain 



We now apply essentially the same procedure to the linear chain, analyzing decoherence 
and deriving equations of motion for the coarse-grained position averages Xj. 



Using the initial condition (|2.17 ) the integrals over q a can be carried out in (|6.3| ) yielding 
an expression for the decoherence functional of the coarse-grained modes, D^[A'(t), A(t)) 
analogous to (|6.9| ). The magnitude of is 



D (d) 



A'(t),A(t) 



exp — 



tf 



, df r dt (flfjKfV, t)&) 

2 Jo Jo K 



(6.20) 



where 



m = A\t) - A(t) 
1 



A(t) 



A'(t) + A{t) 



(6.21a) 
(6.21b) 



The all-important kernel Kj (t',t) is given by 



Kf\t',t) 



N 2 p 2 k B T. 

Q < 

Ah 2 



(tt 2 L I - Ct 2 Q ) M^fT 1 (cos[Qi\ cos[Ot' 



+ sin[fit] sin[Ot'] jCL^M^ 2 (n 2 L I- nj) c 



(6.22) 
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which is proportional to the correlation function (|3.26|) of the classical noise derived from the 
Lagrangian. Thus the "off-diagonal" elements of the decoherence functional decay exponen- 
tially with increasing difference between A'(r) and A{r). Decoherence of coarse-grainings 
with suitably large regions A* at suitably spaced intervals of time is thus achieved. We 
shall return to a detailed discussion of decoherence times as a function of d below, but first 
we consider the equations of motion. 

The probability of a given coarse-grained history a is given by 



p(a)= / 5 AD A{t),A(t) (6.23) 
which in analogy to ( 6.17|) can be expressed as 



p(a) = j 6 A [det (kJ ) /4tt] 5 

x exp{-^ jf*' dt' j tf dt (E\t\A(T)}K.f ym \t\t)E\t,A(T)})) . (6.24) 



S'(t, A(t)] is a linear functional of A{r) analogous to (16.13 ), whose exact form we return to 
below; K^ inv (t',t) is the inverse of K^(t' ,t) defined in flQ2l) . 

Eq ( |6.24| ) shows that the probability of histories is sharply peaked about those obeying 
the equations of motion 

£'(t, A(t)\ = , (6.25) 

but with Gaussian noise causing deviations from this predictability related to the size of 
the kernel K^ mv (t',t). The form of the equations of motion (|6.25| ) could be worked out 
by following the procedure used in the simple example above. However, it is more direct 
to note that the result of that calculation is exactly the equation of motion that would be 
obtained by a completely classical analysis of the coarse-grained dynamics, just as in Q6.17 ) 



for the simple model. The equations of motion £' (t, X (r)] = are therefore equivalent to the 
classical equations derived in Section III and IV. This is an important simplification because 
a straightforward extension of the above analysis would lead to nonlocal terms such as occur 
in ( |6.13D , as we saw in section IIIC In the end, these all cancel to give the simple equations 



of motion exhibited in Section III. This is checked explicitly in Appendix C. 

The size of the kernel (t, t') in (|6.20| ) controls the efficacy of decoherence — the larger 



the kernel the shorter the decoherence time scale. We study it as a function of d holding 
all other parameters fixed, including the temperature of bath T. The kernel is necessarily 
positive |J. A simple measure of its size is the time-averaged trace 

KAd)= lim - [ tf dtTr(K^(t,t)) . (6.26) 

tf-Kx tf JO V ' 

Figure 4 shows JCj plotted as a function of d. Since the form of the kernel is closely 
related to the classical noise correlation function M[A/'(i)A/'(t')] from equation ( |3.26| ), it 
is not surprising that the shape of the result is virtually identical to Figure 3. While Figure 
3 was computed from the form (|3.25|) for the noise, which is different, the difference amounts 
to a constant factor of N 2 and a ^-dependent factor very close to 1. From the two graphs, 
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L = 100 
L = 65 
L = 30 



1000 



10000 



FIG. 4. The time-averaged trace )Cj(d) of the decoherence kernel K.^f\t,t) in units of 



Nk B Tnu 2 /4fr 



vs. 



the coarse-graining d, for Ai = 630 and three L values L = 30, 65, 100. 



The form of these curves is virtually identical to the classical noise correlation function, plotted in 
Figure 3. Note that an overall factor of iV 2 has been divided out to make comparison to Figure 3 
more exact. 
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where the N dependence has been absorbed into the choice of units, it is evident that the 
dependence on d of the two forms of noise is almost exactly the same. 

Since the kernel K/(cf) is a factor of N 2 larger than the classical noise correlation function, 
increasing the level of coarse-graining N exponentially improves decoherence, but actually 
reduces the classical noise [cf. (|5.6|) 1- Because of this, for realistically large values of N, 
the decoherence rate is rapid compared to dynamical time scales even in the completely 
localized case, where the absolute strength of the noise is weakest. One can understand this 
as increasing the absolute level of decoherence, but simultaneously increasing the inertia of 
the coarse-grained variables to resist the increased noise, since the inertia of a group just 
goes like its total mass Nfi. We discuss this trade-off more thoroughly below. 

VII. CONCLUSIONS CLASSICALITY 

For a system with many degrees of freedom, useful dynamical predictions concern reg- 
ularities emerging from coarse-grained descriptions. To be sure, at a fine-grained level the 
system will display the regularities arising from its fundamental equations of motion if it is 
classical, or from the Schrodinger equation if quantum mechanical. However, these regular- 
ities are usually impossible to extract or apply if the number of degrees of freedom is very 
large. The useful predictions arise from much smaller numbers of coarse-grained variables 
correlated in time according to phenomenological equations of motion. 

There are arbitrarily many sets of alternative coarse-grained histories that can be as- 
signed probabilities on the basis of a closed system's initial condition and fundamental 
dynamics. Which of these will exhibit useful regularities in time governed by phenomeno- 
logical equations of motion? How much, and what, coarse-graining is needed to obtain 
useful predictability? This, very roughly, is the problem of characterizing classicality that 
we mentioned in the Introduction. In this conclusion we describe how our results for the 
linear harmonic chain bear on this question. 

Which coarse grained descriptions are predictable is an important problem even in clas- 
sical physics. However, it is especially important in quantum mechanics where probabilities 
can be assigned only to decoherent sets of histories, and two such sets are generally mutually 
incompatible. Further, in a loose sense, the number of sets of coarse-grained histories is 
much larger in quantum mechanics than it is in classical physics. It is therefore important 
to explain within quantum mechanics (and a theory of the cosmological initial condition) 
why we find it useful to employ only a narrow class of the possible coarse-grainings by 
which the universe could in principle be described. Many see this as the central problem in 
understanding quantum mechanics. 

Characterizing classicality involves comparing the utility for prediction of various sets 
of coarse-grained alternative histories. In this paper we have considered only a family of 
coarse-grainings of the simplest linear system exhibiting a continuum description of the kind 
usually found in classical physics — the linear harmonic chain. 

Four parameters characterize the family of coarse-grainings we have considered. The 
atoms of the chain are divided into groups of N, consisting of equally spaced clumps of 
d neighboring atoms. (See Figure 2.) The center of mass coordinates of each group are 
coarse-grained by exhaustive sets of equal ranges A spaced by equal time intervals At. The 
four parameters are therefore N, d, A and At. 
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The probabilities of the individual members of a decoherent set of histories follow from 
the initial state, so that any comparison of the predictability of different coarse-grainings will 
depend crucially on the nature of this initial condition. We assumed an initial condition in 
which the short wavelength modes of the string were in thermal equilibrium at a temperature 
T, while the long wavelength modes were excited well above the level of thermal fluctuations. 
This is an initial state of local but not global thermal equilibrium. 

Both classically and quantum mechanically the probabilities of these decoherent sets of 
alternative coarse-grained histories can be characterized as arising from equations of motion 
augmented by noise. The dynamical time scale of the equations of motion ( 3.12|) is given 
roughly by 

<*» ~ z (t) • (7 ' 1) 

The local coarse-grainings d ~ N have the longest dynamical time scales because they 
correspond to the approximately conserved center of mass momenta of the groups. 

The utility of these equations of motion for prediction depends on the size of the devia- 
tions from the regularity in time they describe that is caused by the noise. The noise forces 
were estimated roughly in eq. ( |5.6| ). (The force is Nfi times the force per unit mass Afjf*). 
For large N, d comparable to N, and L small compared to M.: 

F nmse ~ {k B Tuj^) l l\N/d) l l\L/M) . (7.2) 

With these basic estimates in hand, we can compare the different members of this fam- 
ily of coarse-grainings with respect to three properties bearing on classicality: the rate of 
decoherence, the deviations from predictability caused by the noise, and the computational 
complexity required to use the equations of motion to make predictions. The results are as 
follows: 

• Decoherence: Decoherence and noise are connected. The kernel which governs the 
size of the imaginary part of the influence phase and effects decoherence is the correlation 
function of the noise force [cf. ( |6.15|) and (|6.19|) 1. From ( |6.20|) we can obtain the following 
rough estimate for the decoherence time tdecoh which At must exceed if the coarse-grained 
set of decoherent histories is to decohere: 

tdecoh ~ h/ {F noise A) . (7.3) 

This time increases with the locality of the coarse-graining d because F noise decreases with 
d, as ( |7.2j ) shows. Since the more local coarse-grainings are dominated by low frequency 
modes, the dynamical timescale also increases with d, so that the ratio of the two times 
favors the more local coarse-grainings; combining ( |7.1| ) and (|7.2|) we have 

tdecoh ~ (^£B\(NA\-V* 



tdyn A 



){Nd)-V* , (7.4) 



where Xdb is the thermal De Broglie wavelength introduced by Zurek [£1]] 

Xdb = h/(k B Tfiy/ 2 . (7.5) 
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The decoherence time scale must be less than the dynamical time scale to use the equations 
of motion at all. This can be achieved by taking N or A or both sufficiently large, and 
for both local and non-local coarse-grainings is not a very constraining condition given 



"realistic" parameter values, as was stressed by Zurek J2l|] in simpler cases. For d ~ 1, 
N ~ 10 6 , n ~ WAMU, and T ~ 300°, 

Uecoh/tdyn ~ (l(T 13 cm/A) . (7.6) 

For such coarse-grainings, decoherence is not a major restriction on predictability. 

• Noise: As ( [7.2|) shows, the noise force decreases as d increases for fixed N, that is, it 
decreases as the coarse-graining becomes more local. However, it is not the absolute scale 
of the noise that is important for predictability, but rather its size relative to the dynamical 
force terms F^ yn that occur in the equation of motion. If £ is the characteristic size of the 
excitations of the chain that occur in the coarse graining, then roughly 

F dyn ~ NfiC/t% n . (7.7) 

The size of £ is determined by the initial condition and varies with both iV and d. In 
( |2.16j ) and ( |2.17| ) we assumed an initial condition in which the fine-grained modes above 



a mode number tc were thermally excited, while modes below £c were much more highly 
excited. The connection between fine- and coarse-grained modes given by ( ^.7[ ) and ( [2.9| ) 
shows that the fine-grained modes below lc contribute to coarse-grained modes only when 
d > N(L/£c). Thus, if d lies much below N, we have 



F 

- 1 noise -i 7 7t r /r? n\ 

— 1, d«N, (7.8) 

and the regularities of the equation of motion will be swamped by the noise. 
By contrast, when d ~ iV 

^!££ „ ^VN, d~N (7.9) 

"dyn 

where £t is the characteristic scale of thermal excitations of the mode: 

£ T ~ (ksT/NfiQl) 1 / 2 . (7.10) 

If the size of excitations of long wavelength modes is much greater than that of thermal 
fluctuations, the effect of the noise on the equation of motion will be negligible. Thus 
the local coarse-grainings exhibit more regularity in time if the initial condition has this 
property. For the realistic string described above with iV ~ 10 6 , typical thermal excitations 
of the L = 10 mode would be £t ~ 10~7 m. 

• Computational Complexity: Even when coarse-grainings decohere so that probabilities 
can be assigned to histories, even when noise is negligible so that histories exhibit the 
regularities in time summarized by classical equations of motion, coarse-grainings can be 
distinguished by the effort needed to calculate these regularities. The number of operations 
Ns necessary to evolve the equations of motion over a time interval T is roughly proportional 
to 
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N s oc (M/t dyn )T ~ (LujT /d) 



(7.11) 



Thus prediction becomes easier as both N and d increase. For "realistic" coarse-grainings 
where N ~ 10 6 there can be a significant difference between the local coarse-grainings with 
d ~ 10 6 and the nonlocal ones with d ~ few. 

In summary, decoherence, resistance to noise and computational simplicity all favor local 
coarse-grainings over nonlocal ones — all these comparisons being contingent on large N and 
an initial condition in which long wavelength modes are more excited than short wavelength 
ones. 

These quantitative results for the linear harmonic chain support the heuristic arguments 
for the predictability of more general kinds of quasiclassical variables that were sketched in 
the Introduction. In our family of coarse-grainings, the ones more useful for prediction are 
the more local ones associated with an approximately conserved quantity. Our analysis of 
the harmonic chain is a step towards a more realistic analysis of classicality in at least three 
ways: (1) It considered a system which permits a continuum approximation of the kind 
usually found in classical physics. (2) It employed a system-environment split which follows 
from the coarse-grainings needed to realize that approximation rather than being posited 
ad hoc in terms of fundamental coordinates. (3) Different coarse-grainings were compared 
quantitatively with respect to decoherence, noise and computational complexity. 

However, these positive features should not obscure how short this analysis falls from 
the kind of treatment of classicality envisaged by We considered only linear 

interactions, not a realistic Hamiltonian. We did not compare all possible sets of alternative 
coarse-grained histories, but only a four-parameter family of them. We did not propose a 
unified quantitative measure for classicality, but rather dealt separately with some of its 
attributes: decoherence, resistance to noise, and computational simplicity. We did not start 
from the initial condition of the universe nor exhibit the important role played by gravity in 
creating the conditions for local equilibrium while ensuring the absence of global equilibrium. 
Rather we assumed these properties in our initial condition. Future analyses will do better. 
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APPENDIX A: SYSTEM-ENVIRONMENT SPLITS 

In this appendix we describe the circumstances under which — for a given coarse-grained 
set of alternatives — the Hilbert space 7i of a closed system can be written as a tensor 
product 7i s ® H e where Ti s contains the quantities followed by the coarse-graining and TL e 
contains the quantities that are ignored. Such a tensor product factorization is called a 
system-environment split. 
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We begin by considering a Hilbert space 7i and a set of alternatives at a single moment of 
time represented by an exhaustive set of orthogonal projection operators {P a }, a = 1, 2, • • • 
satisfying 

P a Pp = 5 aP Pp, y £P a = I. (Al) 



We seek to write 



such that 



n = n s ®n e (A2) 



P a = P s a ®I\ a = l,2,.-.. (A3) 

The decomposition is then a system-environment split. If 7i s were to contain just the 
quantities followed by the coarse-graining we would naturally impose dim (P*) = 1, where 
dim (P) is the dimension of the subspace projected on by P. However, other notions of an 
environment can be useful in which (dimP*) > 1, and we shall consider the general case. 

Only two simple mathematical facts are needed to analyze the above question. First, the 
decomposition (|A3|) requires the relation between dimensions 



dim (P Q ) = dim (P*) dim (H e ) , (A4) 

and, clS db special case 

dim (H) = dim (H s ) dim (H e ) . (A5) 

Second, two separable Hilbert spaces of a given dimension are isomorphic — a consequence 
of the fact that they both have countable bases. 

Realistic cases have dim (TC) = oo, but we pause to note some evident results from ( |A4| ) 



when dim is finite. Then all the dim (P a ) are finite. A system-environment split is not 
always possible, only when the dim (P a ) are all divisible by a common factor. In particular, 
if dim (P*) is required to be unity, then a system-environment split is possible only when all 
the P's have the same dimension. 

When dim (7i) is infinite, there are a number of subcases which are convenient to treat 
separately. The most important of these is when dim (P a ) = oo for all a. The requirements 
( |AT| ) imply that the P a all commute and can be simultaneously diagonalized. Let {\i)}, 
i — 1, 2, • • • be a basis in which they are all diagonal. Then to each P a there is a subset of 
these basis vectors spanning the corresponding subspace. We write 

P a = £|z)(z|. (A6) 

However, it is then a simple matter of relabeling to define an isomorphism between the 
infinite dimensional Hilbert spaces 7i and 7i s ® 7i e . We write 

\i) = \a,A) (A7) 
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where i ranges over the infinity of states in P a and A = 1, 2, • • • is another labeling of them. 
This relabeling defines the tensor product 7i s <8> 7~C e - Operators acting only on 7i s have the 
form 

(a'A'\0\aA) = (a'\O s \a) 8 AA > (A8a) 
while those on H. e have matrix elements 

(a'A'\0\aA) = 5 a , a (A'\O e \A) . (A8b) 

In particular, 

P a = P s a ®P = (\a)(a\)®I e . (A9) 

and dim (P*) = 1. We have constructed a system-environment split defined by the coarse- 
graining {P a }. 

If the condition dim (P*) = 1 is relaxed it is possible to define many other system- 
environment splits for this coarse-graining. One simply relabels including more states in 7i s , 
viz. 

\i) = \a,A) (A10) 

such that the subspace P^ contains several different values of a. 

Many calculations use a system-environment split of this kind. For example, in studies 
of Brownian motion the labels a correspond to the coordinates of the massive particle and 
A to the coordinates of the particles of the bath. The coordinates of the bath are ignored, 
but typically the coordinates of the Brownian particle are followed only to some accuracy. 
Thus, for a given choice of ranges, coordinates other than those in the bath are ignored 
corresponding to dim(P^) > 1. There is no unique system-environment split. 

The key to the above construction is that the relation ( |A4| ) is easily satisfied because 
dim (P a ) and dim (7i) are both infinite. Finite dimensional members of the set of alternatives 
are obstacles to a system-environment split. If dim (P a ) is finite for some a, then can 
only be satisfied if dim (P*) and dim (7i e ) are both finite — already a restrictive condition. 
Furthermore, the dimensions of all the finite dimensional P a must be divisible by common 
factor — dim (7i 6 ) — and this is not always possible. The only case in which the condition 
dim (P*) = 1 could be enforced is if all the finite dimensional P a have the same dimension. 
Clearly, a system-environment split is generally possible only when the dimensions of all of 
the Pq, are infinite. 



APPENDIX B: MATRICES FOR THE CHAIN OF OSCILLATORS 

In this appendix we complete a calculation begun in Section IIIC This is the explicit 
demonstration that when the variables in the action are changed according to ( p.9| ) and 
fl2.14| ), the result yields the same equations of motion as are obtained by making these 
changes in the fine-grained equations of motion directly. Specifically we check that when 
( p.7| ) is used to eliminate the g's from ( |3.19|) it yields ( p.12 ). 
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Since the coarse-grained modes A^\t) are all decoupled from each other, and interact 
with separate collections of high-frequency modes, we can consider them one at a time, 
which somewhat simplifies our notation. In this case, the matrix S reduces to a diagonal 
matrix with elements making up a single vector with d components S a and T is a d by d — 1 
matrix with components T ab , where a ranges from to d — 1 and b ranges from 1 to d — 1. 
Thus M ss and V ss become scalars, and M ST , M TS , Y ST and V T5 are vectors in the space 
of the g's. These matrices have the simple form 

S = —, S a = 0, for a > (Bl) 

c 

and 

T b = ~— , T ab = S ab , for a > (B2) 
c 



where the numbers c a are the coefficients c \j, defined by (|2 .11 ) and ( |3.6| ), forming the com 



ponents of the vector c The matrix u is diagonal with components u a and the mass is a 
constant /i. 

With these definitions, the matrices defined in (|3.17a| - |3.17b| ) are 



M : 



ss 



I c 

(M TT ) W = y,5w + ^ 



2 : 

C b C*y 



c 



2 ' 



(M S T)b = -Vr% = (M TS )b, 
I c o| 



V S s = 



Cb 



(V 3T ) b = -^oVW = (VTs)b. (B3) 



Co 



In these expressions when the matrices M§s, etc., have been reduced to scalars or vectors 
we have changed the notation in what we hope is an obvious way. With (|B3|) we can now 
show that equation (|3.19|) is identical to equation ( |3.12| ) aside from a multiplicative factor. 
With a little algebra we find that 



(M ss - M ST M^M TS 



|c | 2 + Eb \cb\ 2 , 



(V ss - MstM^Vts) = ( ,„ , 2 JL | | 2 ^1 



-(V S T - M ST M T l T V TT ) b = ^ 



|c | 2 + T,b \Cb\ 

/' 



|co| 2 + J2b' |Cfe'| 2 / 

Since uq = Ql, plugging these back into ( |3.19| ) simply yields the equation 



CbH-ujt). (B4) 



■ l2 JL , |2 [Af{t) +O|4 d) (0) = | , 2 / v | J A/(t), (B5) 
co r + Eft C6 \l c ol + E& Q, 2 / 
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that is, simply equation fl3.12|) multiplied by a constant factor. This factor is independent 
of the choice of coarse-grained mode L or the coarse-graining d: 



fi \ / /i 



|co| +T,b\cb\ 2 J \\c \ 2 + cc 



N/i . (B6) 



APPENDIX C: TRANSFORMING THE EQUATION OF MOTION 

In Section III we derived the classical equation of motion for the coarse grained variables 
Ai(t) utilizing two different ways of eliminating the environmental coordinates q a . First, 
using the fact that the g's were themselves fine-grained coordinates, we solved the classical 
equations of motion for the g's and substituted the solution into the fine-grained equations 
of motion. Second, we derived the equations of motion from the Lagrangian written in terms 
of the A's and g's and then solved the for the g's to eliminate them. In this appendix we 
complete the demonstration that the ostensibly different equations for the A's that result 
are, in fact, equivalent. 

Suppose we have an equation of motion of the form (|3.9|) for a single coarse-grained mode 
Al, which in this appendix we call Sl = 0. We wish to write it in a new form (|3.22j ) which 
we call S' L = by a transformation of the form: 

S' L {t,A L {r)\ = C (s L {t,A L {r)\ + jf* dt'G L (t,t')S L (t' \A l (t)^ = Af L (t), (CI) 

with, of course, a transformed noise function 

Af' L (t) = C (Af L (t) + jf* dt'G L (t,t')Af L (t'f) . (C2) 

Here, C is a positive constant and is a particular Green's function for the equation 

Sl = 0, that is in operator shorthand SlGl = I ■ In a similar shorthand we refer to the 
transformation ( |CT| ) as C(I + G). Clearly, the solutions to the two equations ( ^.9[ ) and ( C 1| ) 
will be the same only if (/ + G) is invertible. This will be true if the equation 

h(t)+ f ' G L (t,t')f L (t')dt' = (C3) 
Jo 



can only be solved by fhif) — 0. From ( |C3| ) it is clear that any solution fi(t) must have 
/l(0) = 0, and is a solution to a second-order linear equation. 



In the case of the chain of oscillators, the original form of the equation is (3.12) and the 



transformed equation is ( [3.22j ). From Appendix B we see that the positive constant is 



C = , /, J = Nim, (C4) 



Co + cc 



and from (ET22Q and (M§ the kernel G L (t,t') is 



G L {t, t') = -^-c(Q 2 L I - fij) M^fi" 1 sin[n(t - t')]M^ /2 c, (C5) 

Co 
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where Q 2 L is the constant uo 2 and fig is the diagonal matrix uo restricted to the g's, with 
diagonal elements u 2 . 

From ( Q5| ) we know that Gi{t,t) = for all t, which implies 



= is™*''*' (C6) 

and hence (dfL/dt)(0) = 0. Thus, the only solution to ( |C3"D is fhif) = and therefore (I+G) 
is indeed invertible. 

Now we need to show that C{I + G)Afi(t) = Af' L (t) for our harmonic chain. The two 
forms of the noise are given by (|3.23a|) and (|3.23b| ), respectively. Let ft 2 be the transformed 
frequency matrix defined by ( |3.21| ), with eigenvalues z/f and orthonormal eigenvectors Vk, 
and On be the diagonal frequency matrix with eigenvalues u 2 . Then the requirement (|C2|) 
can be written 



Nfic{Q 2 L I - H|) (cos(Ogt) Aq(0) + Qq 1 sin(n Q t)A'q{0) 

+ p^SM^ /2 (Q 2 rI - nDQ- 1 /"*sin(fi(t - £'))M^ /2 c 
I c 1 - •/() 

xc{Q 2 L I - n 2 Q ) (cos(O Q t)Ag(0) + Oq 1 sin(O Q t)A^(0)) dt' 
= Nfic(n 2 L l - n 2 ) (cos(rit)M^Ag(O) + rr 1 sin(ftt)M^A<f(0)) . (C7) 

By writing the vectors in terms of eigenvectors of f2 2 we can break the integral in ( p7| ) into 
a sum over many integrals having the forms 



/ sin( u (t — t)) sin(ut) dt — — - [u sm(ujt) — lu sm(ut)] , (C8) 

Jo v z — UJ Z 



and 



sin(z/(t - t')) cos(ut)dt = - - - [cos(cjt) - cos(z4)l . (C9) 

V 1 — ijj 1 



Substituting these values into (|C7|), the result we wish to show becomes 
5> (Ql - u 2 a ) fcos(cj a t)Ag a (0) + — sin(w a t)Ag a (0) 



+ I ^JllcM-TT V k \ 2 C a - 



_ cos(w a t)Ag a (0) 

- cos(z/ fc t)Ag a (0) + — sin(^ a t)Ag a (0) - — sin(u k t)Aq a (0) 
= E(cM^ /2 *4)(^i - v\) fcos(z/ fc t)(4M^Ag(0)) + - sin^X^M^A^O))) . (CIO) 

We wish this to hold at arbitrary times t, which implies that each frequency must be equated 
separately. This requires that the following two conditions hold true: 

I + -^cM^ /2 (fi 2 I - n 2 )(Q 2 - u 2 a I)- 1 M~ 1 T /2 c = , (Clla) 
I I 
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(4M^Ag(0)) + J-(^M-y 2 c)^c a ^f ^ Ag a (0) . (Cllb) 

l C 0| a l^fe _ W J 

The first condition ( |Cllaj) must hold for all a. The last condition (|Cllb|) must hold for 
arbitrary initial vector Ag(0) and all k. 

We can readily evaluate (|Clla ) from the matrix definitions ( p3|) . The inverse matrix is 

(fi 2 - u 2 a I)- 1 = M^(V TT - ulM TT )- 1 'M^}, (C12) 

and the matrix (Vtt — w^M-jt) has the simple form My = fcSij + g^g^ and hence can be 
analytically inverted (even though one of the fi vanishes). Carrying out this inverse and 
performing the sum, we see that the a dependence of QUI lap does indeed drop out, and the 
equation is satisfied. 

The second condition ( Cllb ) looks even more difficult to evaluate, since we have no 
explicit expressions for the eigenvalues v\ and eigenvectors v^. However, these are not 
necessary. Note that if v*. is an eigenvector of f2 2 = M^/ 2 \ VtM^ 2 with eigenvalue z/f, 

1 /2 — * 1 

then ~M.rp T Vk is an eigenvector of V^rM^ with the same eigenvalue. Using the definitions 
), this implies that 



Since c is an eigenvector of Mtt, we can readily show 

(vkM^c) = - -. ^ o| ' (v k Mjffc) . (C14) 



By combining ( |C13| ) and ( |C14j) , the condition ( |Cllb| ) is immediately proven. Hence, the 
transformed classical noise is the same as the noise derived from a straightforward Lagrangian 
calculation. 

This result has immediate implications for the quantum probabilities ( |6.24 ). The inte- 
grand in the exponent of that expression is proportional to 

£' L (t, A L (r)]M[Af (t) Af (t')] inv £' L (t' , A L {r)] (C15) 

summed over all L. This is all in terms of the equation of motion form (|3.22|) and noise 
correlation function But since E' L = C(I + G)S L and M[Af'(t)Af'(t')] inv = C' 1 ^ + 

G)- 1 M[A/(t)A/(t')] inv C~ 1 (J + G)~\ we see that the transformation C(I + G) cancels out, 
and 

8' L {t,A L {r)}M[Af\t)Af\t')r£' L {t\A L {r)} = S L (t,A L (r)}M[Af(t)Af(t')rS L (t',A L (r)] , 

(C16) 

i.e., even in the quantum case one can use the simpler form of the classical equation of 
motion and noise. 
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